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Supersymmetry has been often invoked as the new physics that might reconcile the experimental 
muon magnetic anomaly, a M , with the theoretical prediction (basing the computation of the hadronic 
contribution on e + e~ data). However, in the context of the CMSSM, the required supersymmetric 
contributions (which grow with decreasing supersymmetric masses) are in potential tension with 
a possibly large Higgs mass (which requires large stop masses). In the limit of very large mj, 
supersymmetry gets decoupled, and the CMSSM must show the same discrepancy as the SM with 
a M . But it is much less clear for which size of mj does the tension start to be unbearable. In this 
paper, we quantify this tension with the help of Bayesian techniques. We find that for mh > 125 
GeV the maximum level of discrepancy given current data (~ 3.2 a) is already achieved. Requiring 
less than 3 a discrepancy, implies mi < 120 GeV. For a larger Higgs mass we should give up either 
the CMSSM model or the computation of based on e + e~; or accept living with such inconsistency. 



I. INTRODUCTION 

The magnetic anomaly of the muon, a M — \{g~ 2) M has 
been a classical and powerful test for new physics. As it 
is known, the present experimental value and some of the 
theoretical determinations of o M show a remarkable dis- 
crepancy, suggesting physics beyond the Standard Model 
(SM) to account it. However, the situation is still uncer- 
tain, due essentially to inconsistencies between alterna- 
tive determinations of the contribution coming from the 
hadronic vacuum-polarization diagram, say c^ad^- 

This contribution can be expressed in terms of the to- 
tal hadronic cross section e + e _ — > had. Using direct 
experimental data for the latter, one obtains a final re- 
sult for a^ t , which is at more than 3 a from the current 
experimental determination [lj, namely 

<5a M = a c ^-af l = 25.5 ± 8.0 x 1(T 10 (1) 

(the quoted error bars are 1 a). This discrepancy has 
been often claimed as a signal of new physics. Obviously, 
if one accepts this point of view, the discrepancy should 
be cured by contributions from physics beyond the SM. 

Admittedly, such claims are too strong. We are quite 
aware of past experimental observables in apparent dis- 
agreement with the SM prediction, which have eventually 
converged with it. This has occured due to both experi- 
mental and theoretical subtleties that sometimes had not 
been fully understood or taken into account. As a mat- 
ter of fact, the experimental e + e~ — > had cross section 



* E-mail: maria.cabrera@uam.es 
t E-mail: alberto.casas@uam.es 

tE-mail: rruiz@ific.uv.es 

§E-mail: r.trotta@imperial.ac.uk 



exhibits some inconsistencies between different groups of 
experimental data. Using only BABAR data the discrep- 
ancy reduces to 2.4 a, while without it the discrepancy 
becomes 3.7 a, [Tj. The inconsistency is specially noto- 
rious if one considers hadronic r decay data, which are 
theoretically related to the e + e~ — > had cross section. 
Using just r-data the disagreement becomes 1.9 a, [T], 
[2J. Although the more direct e+e~ data are usually pre- 
ferred to evaluate a^ M , these inconsistencies are warning 
us to be cautious about the actual uncertainties involved 
in the determination of a^ M . 

If one takes the discrepancy between theory and exper- 
iment shown in eq.Q as a working hypothesis, one has to 
consider possible candidates of new physics able to pro- 
vide the missing contribution to reproduce a° xp . The 
Minimal Supersymmetric Standard Model (MSSM) is 
then a natural option. We will consider here the simplest 
and most extensively analyzed version of the MSSM, 
namely the so-called constrained MSSM (CMSSM), in 
which the soft parameters are assumed universal at a high 
scale (Mx), where the supersymmetry (SUSY) breaking 
is transmitted to the observable sector, as happens e.g. 
in the gravity-mediated SUSY breaking scenario. Hence, 
our parameter space is defined by the following parame- 
ters: 

{0} = {m,M,A,B,fx,s} . (2) 

Here m, M and A are the universal scalar mass, gaug- 
ino mass and trilinear scalar coupling; B is the bilinear 
scalar coupling; [i is the usual Higgs mass term in the su- 
perpotential; and s stands for the SM-like parameters of 
the MSSM, i.e. essentially gauge and Yukawa couplings. 
All these initial parameters are understood to be denned 
at M x - 

The main supersymmetric (CMSSM) contributions to 
a M come from 1-loop diagrams with chargino-sneutrino 
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and neutralino-smuon exchange [3J- In general, these 
contributions, say 5 MSSM a M , are larger for smaller super- 
symmetric masses and can be just of the right magnitude 
to reconcile theory and experiment (thus constraining the 
CMSSM parameter space). 

In section [TT] we show the potential tension between 
the requirement of suitable SUSY contributions to the 
muon anomaly and a possibly large Higgs mass. In sec- 



tion III we quantify such tension as a function of m/,, 
with the help of Bayesian techniques. In section |IV| we 



show how the probability distributions of the most rel- 
evant parameters (universal scalar and gaugino masses, 
and tan (3) change with increasing rrih . Finally, in section 
5 we present our conclusions. 



II. HIGGS MASS VS. G-2 

It is well known that in the MSSM the tree-level Higgs 
mass is bounded from above by Mz, so radiative correc- 
tions (which grow logarithmically with the stop masses) 
are needed to reconcile the theoretical predictions with 
the present experimental lower bound, rrih > 114.4 GeV 
(SM-like Higgs). Roughly speaking, a Higgs mass above 
130 GeV requires supersymmetric masses above 1 TeV. 
In this regime one can expect SUSY to be decoupled, so 
that the prediction for <x M becomes close to a^ M . Hence, a 
large Higgs mass in the MSSM would necessarily amounts 
to a > 3 a discrepancy between the experimental and 
the theoretical values of a M (evaluated via e + e~ — > had). 

The main goal of this paper is to quantify the tension 
between rrih and in the context of the CMSSM. This 
is useful since it allows to put an educated upper bound 
on the Higgs mass, which will depend on the discrepancy 
one is ready to tolerate. Conversely, it tells us from which 
minimum value of rn c h * p we will have to give up either the 
CMSSM assumption or the theoretical evaluation of 
via e + e~ — > had (with the quoted uncertainties). 

For the sake of the discussion, we will give now some 
approximate analytical expressions for m/, and <5aJf SSM . 
In the MSSM the tree-level squared Higgs mass plus the 
one-loop leading logarithmic contribution is given by 
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Here tan/3 is the ratio of the expectation values of the 
two MSSM Higss fields, tan/3 = (H u )/(H d ); m t is the 
(running) top mass and m t ~ is the geometrical average of 
the stop masses. Besides, 



X t = A t + ii cot , 



(4) 



where A t is the top trilinear scalar coupling, and M| is 
the arithmetical average of the squared stop masses. All 
the quantities in eqs.([3]), Q are understood at low en- 
ergy (for more details see e.g. ref . |3HTU] ) . Subdominant 



terms not written in eq.Q can be important for a pre- 
cise determination of rrih, and we have included them in 
the numerical analysis. The previous equations tell us 
how to/j grows with increasing supersymmetric masses 
and also with increasing tan/3. Besides, the contribution 
associated to the stop mixing (second term within the 
square brackets in eq.([3])) is maximal at X t — y/6Ms- 

On the other hand, as mentioned above, the super- 
symmetric contribution to the muon anomaly, <5 SUSY a At , 
arises mainly from 1-loop diagrams with chargino- 
sneutrino and neutralino-smuon exchange. This contri- 
bution increases with increasing tan /3 and decreasing su- 
persymmetric masses. See refs. |lH - [T5] . 

Although the analytical expressions are complicated, 
one can get an intuitive idea of the parametric depen- 
dence by considering the extreme case where the masses 
of all supersymmetric particles are degenerate at low en- 
ergy 1 : Mi = M 2 = n = mjxL = mp,R = m 5 = A/susy • 
Then M, 
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g\ tan /3 sign(M 2 //). 



(5) 



Examining the approximate expressions ([3| and ([s]), it is 
clear that a large rrih and a large 5 SUSY a M will be more 
easily obtainable (and thus compatible) for larger tan /3. 
On the contrary, the larger the supersymmetric masses 
the larger rrih but the smaller 5 SUSY a M , and this is the 
origin of the potential tension. 

However, it is difficult from the previous expressions 
(or the more sophisticated ones) to conclude for which 
size of rrih does the tension start to be unbearable. The 
reason is that a particular value of the Higgs mass, say 
rrih = 120 GeV, can be achieved through eq.([3]) with 
different combinations of tan/3, stop masses and X t . Be- 
sides, there are many ways, i.e. very different regions 
in the MSSM parameter space, in which these quantities 
can have similar low-energy values. Still, the correspond- 
ing contribution 8 can change significatively from 
one region to another. Unless one performs a complete 
scan of the parameter space one cannot conclude that the 
required value of 5 SUSY a M is unattainable for rrih = 120 
GeV. On the other hand, if it is attainable, but only in 
an extremely tiny portion of the parameter space, this 
implies a tension between the two observables since the 
consistency between rrih and requires a severe fine- 
tuning. And it is possible, in principle, to quantify such 
tension. 

In the analysis we have included two-loop leading cor- 
rections for the Higgs sector [TTHSTj . <5 SUSY a /J was com- 
puted at full one-loop level adding the logarithmic piece 
of the quantum electro-dynamics two-loop calculation 
plus two-loop contributions from both stop-Higgs and 
chargino-stop/sbottom [T3]- The effective two- loop effect 



1 This limit is often used because of the simplification of the for- 
mulae it implies. However, it is unachievable in the CMSSM. 



3 



due to a shift in the muon Yukawa coupling proportional 
to tan 2 j3 has been added as well [14] . 

Next we expound how a systematic analysis of this 
kind can be done with the help of Bayesian techniques. 
This will allow us to quantify the tension between rrih 
and o„ as a function of to/,. 



III. QUANTIFYING THE TENSION BETWEEN 

m h AND o M 

Le us start by recalling some basic notions of Bayesian 
inference. We refer the reader to [S2J [33] for further de- 
tails. For a model defined by a set of parameters 0, the 
posterior probability density function (pdf) of a point in 
parameter space, {9}, given a certain set of data, is de- 
noted by p(0|data) and it is obtained via Bayes theorem 
as 



p(0|data) 



p(data|0) p{9) 
p(data) 



(6) 



Here p(data|0) is the likelihood function (when consid- 
ered as a function of 9 for the observed data) 2 . p(9) is 
the prior, i.e. the probability density that we assign the 
points in the parameter space before seeing the data (in 
the context of Bayesian inference, the prior for a new cy- 
cle of observations can be taken to be the posterior from 
previous experiments). Finally, p(data) is a normaliza- 
tion factor, sometimes called the evidence. It is given 
by 



p(data) = d9 p(data|0) p(9) 



(7) 



i.e. the evidence is the average of the likelihood under 
the prior, and thus it gives the global probability of mea- 
suring the data in the model. 

When two different models (or hypotheses) are used 
to fit the data, the ratio of their evidences gives the rel- 
ative probability of the two models in the light of the 
data (assuming equal prior probability for both). For 
an application to model selection in the context of the 
CMSSM, see [25]. 

In order to quantify the tension between rrih and o^, 
following Ref. [35] we separate the complete set of data 
in two subsets: 



{data} = {9, D}. 



(8) 



Here 9 represents the subset of observations, whose com- 
patibity with the rest of the observations, D, (which are 
assumed to be correct) we want to test. In our case, @ is 



2 Frequentist approaches, which are an alternative to the Bayesian 
framework, are based on the analysis of the likelihood function in 
the parameter space; see ref. |24| for a recent frequentist analysis 
of the MSSM. 



the experimental value of a^, whereas D is given by all 
the standard electroweak observables, B- and D-physics 
observables, limits on supersymmetric masses, etc (for 
the complete list of experimental data used in this pa- 
per, with references, see Table 2 of [37J ). D includes also 
the value of rrih that we are probing, and thus provision- 
al assumed to be the actual one. Hence, we will not 
consider any experimental error in the value of mh, just 
the uncertainty associated to the theoretical calculation 
(estimated as ±2 GeV). Now we construct the quantity 
p($)\D), i.e. the probability of measuring a certain value 
for @, given the known values of the remaining observ- 
ables, D, 



p{@\D) 



p{9, D) 
pjD) ' 



(9) 



Here p{2#, D) = p(data) is the joint evidence, as given 
by Eq. (7]), i.e., the global probability of measuring both 
sets of data at the same time, and p(D) is its equivalent 
but just for the D subset. The latter is a normalization 
factor which will soon cancel out. 

Now, the consistency of %)° hs (the measured muon 
anomaly) with the rest of data, D, in the context of the 
model (CMSSM), can be tested by comparing p{@ ohs \D) 
with the value obtained using different values of in 
particular the one that maximizes such probability, say 
(assuming the same reported error at the new cen- 
tral value). This gives a measure of the likelihood of the 
actual data, ^ obs , under the assumption that the model 
is correct: 



p(@ ohs \D) _ p(3i obs ,D) 
p(3> max \D) ~ p(@ max ,D) 



^(@ obs \D). (10) 



J?(3> ohs \D) is analogous to a likelihood ratio in data 
space, but integrated over all possible values of the pa- 
rameters of the model. Therefore, it can be used as a 
test statistics for the likelihood of the data being tested, 
@ ohs , in the context of the model used (the CMSSM). 
Such test was called Jzf— test in Ref. [26]. Note that, 
as mentioned above, the p(D) factor cancels out in the 
expression of _5f (2# ohs \D), which is simply given by the 
ratio of the joint evidences. 

In our case, the value of 3> max depends on the value 
of TO/i probed. For very large m/,, say rrih > 135 
GeV, SUSY must decouple, so £F max should approach 
the SM prediction. Hence, in this limit one expects 
Jf(^° hs \D) to show a 3.2 a discrepancy; in other words, 
-21nJ^(S> obs |L>) -> 3.2 2 . However, the expression (fTol) 



allows us to evaluate this likelihood for any intermediate 
value of TO/j, and so we can evaluate how quickly this 
limit is reached as a function of the assumed value for 
m h . 

For the numerical calculation we have used the 
MultiNest [28rt30] algorithm as implemented in the 
SuperBayeS code [3TH33] . It is based on the framework 
of Nested Sampling, recently invented by Skilling [541 135] . 
MultiNest has been developed in such a way as to be an 
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FIG. 1: The — 21n^f test statistics, as defined in Eq. (10 1, as 
a function of the assumed value for mi, in the CMSSM frame- 
work with logarithmic (blue, lower dotted line) and flat (vio- 
let, upper dotted line) priors. Horizontal lines denote thresh- 
olds of 2.5 cr, 3 a and 3.2 a discrepancy. 





26 




24 




22 




20 


(3 


18 


►J 




rN 


16 




14 




12 




10 



5 o- 



4tr 



3,r 



Log Priors 
Flat Priors 



110 115 120 125 130 135 140 
m h (GeV) 

FIG. 2: As in Fig. |III| but (artificially) assuming an improved 
experimental determination of a M , so that the SM discrepancy 
becomes 5 cr. Now the horizontal lines denote thresholds of 
3 cr, 4 a and 5 a discrepancy. 



extremely efficient sampler even for likelihood functions 
denned over a parameter space of large dimensionality 
with a very complex structure as it is the case of the 
CMSSM. The main purpose of the Multinest is the com- 
putation of the Bayesian evidence and its uncertainty 
but it produces posterior inferences as a by-product at 
no extra computational cost. 

Fig. 1 shows the value of — 21nJz? (the analogous of 
the usual \ 2 ) 10 r different values of the Higgs mass, 
m h (GeV) = 115,120,125,130,135, and for two differ- 
ent choices of initial priors for the CMSSM parameters, 
namely log prior (red line) and flat prior (blue line) . The 
precise shape of the log and flat priors used here is the 
one derived in ref. |27j . to which the reader is referred, 
which take into account the likelihood associated to the 
electroweak breaking process. The horizontal error bars 
reflect the uncertainty in the theoretical computation of 
to/j in the MSSM, while the vertical error bars come from 
sources of error in the computation of Jzf, mainly the 
numerical accuracy of the evidence returned by Multi- 
Nest. Lines of conventional confidence levels thresholds 
in terms of number of a are shown as well for comparison. 

From the figure we see that the likelihood of the ex- 
perimental value of approaches asymptotically the ex- 
pected 3.2 cr discrepancy for large values of mh, for both 
types of priors. As mentioned above, this is logical and it 
represents a nice cross-check of the reliability of the whole 
procedure. Besides, Fig. 1 tells us how fast this conver- 
gence is reached as rrih increases. And, as a matter of 
fact, the convergence is very fast. At = 125 GeV the 
maximum level of discrepancy is already achieved, indi- 
cating that SUSY has decoupled, and thus the prediction 
for a M coincides with the SM one. If we require less than 
3 a discrepancy, we need to/j < 120 GeV. This is a predic- 



tion of the CMSSM provided we accept the calculation 
of dp based on e + e~ data. For a larger Higgs mass we 
should give up either the CMSSM model or the compu- 
tation of et M based on e + e~; or accept living with such 
inconsistency. These are the main conclusions of this pa- 
per. They stem directly from Fig. 1. Let us also note 
that, even assuming a Higgs mass as low as it can be, the 
minimum level of discrepancy is about 2.5 a. However, 
most of this tension with comes from b — > s, 7 data 
[26], rather from the value of the Higgs mass. This can be 
checked by repeating the analysis excluding all the exper- 
imental information (except Mz and the assumed Higgs 
mass). The resulting plot is similar to that of Fig. 1, 
except the m/, = 115 GeV point, which shows a ~ 1.5 a 
discrepancy. 

It is an interesting exercise to compute how our conclu- 
sions would change if a M became more precisely measured 
in the future (keeping the same central value). If one con- 
tinued to assume the theoretical evaluation of based 
on e + e~ data, the signal for new physics would obvi- 
ously become stronger. In this case, the tension between 
a large Higgs mass and the experimental would get 
more unbearable. We have done this excercise, by chang- 
ing (artificially) the experimental uncertainty of a^ xp , so 
that the discrepancy with the SM result be 5 cr, some- 
thing that could happen in the next years. Now, in the 
context of the CMSSM, the value of -21n J^(^ obs |L>) 
must approach asymptotically such 5 a discrepancy, and 
this is indeed what we observe, as shown in Fig. 2. In 
this hypothetical situation, a Higgs mass above 120 Gev 
would imply a discrepancy larger than 4 a with the muon 
anomaly in the context of the CMSSM. Actually, the 
present lower bound, m/j > 114.4 GeV, would already be 
inconsistent with the muon anomaly at the 3 a level. 

This gives a fair idea of the tensions within the CMSSM 
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to accommodate a value of a M as the measured one (bas- 
ing the theoretical calculation on present e + e~ data). 

IV. PROBABILITY DISTRIBUTIONS FOR 
SUPERSYMMETRIC PARAMETERS 

It is also interesting to investigate the probability dis- 
tributions of the CMSSM parameters for various assumed 
values of the Higgs mass. Figure 3 (upper panels) shows 
the marginalized probability distribution functions (pdfs) 
of m, M assuming a value of rrih = 115, 120, 125 (GeV), 
as well as adding in all present-day constraints mentioned 
above. The location of the peak in the posterior pdf in- 
creases with the assumed Higss mass since, as mentioned 
in section |nj in the MSSM a large rrih requires large ra- 
diative contributions, which grow logarithmically with 
the stops masses. This happens even though large values 
of m and M are penalized both for a natural electroweak 
breaking (see ref. [271 [36]) and by the need of a sizeable 
<5 SUSY a M . The model "prefers" to reproduce m h at the 
cost of not reproducing a M rather than viceversa. Note 
here that for increasing soft masses the discrepancy of 
with the experimental value approaches 3.2 a, but 
if the soft masses are not large enough, the discrepancy 
associated to rrih would be much more severe. 

Fig. 3 (lower panel) shows the pdf of tan /3 for rrih — 
115, 120, 125 (GeV). Its shape is the result two compet- 
ing effects. On the one hand, large values of tan/3 are 
severely penalized for the electroweak breaking [STJ ES] • 
On the other hand, the need of a sizeable 5 SUSY a M favours 
large tan/3 (see the approximate expression ([5])). Fig. 3 
shows the balance between these two effects. [The Higgs 
mass increases also with tan /3, but the effect is only im- 
portant for small values of tan/3, see eq.(|3])]. Now, since 
for larger rrih the soft masses are larger, with the side- 
effect of suppressing one might expect that the 
preferred value of tan /3 increases with m/,, to compensate 
this in eq. ([5]). However, this effect is not very important, 
as it is apparent in Fig. 3. To understand this, let us ap- 
proximate (for the sake of the argument) M susy ~ rrie L R 
in eq. ^ and use [37], [33] 

mj L ~ to 2 +0.54M 2 , (11) 
m\ R ~ to 2 +0.15M 2 , 
mf ~ 3.36A/ 2 + 0.49m 2 — 0.05A 2 — 0.19AM + mf. 

Since rrih increases (logarithmically) with m 2 , while 
<5 SUSY a M is suppressed by m\ , it might seem that the 
most efficient way to reproduce both is to increase M 
rather than m (note the different dependences on M in 
eqs. (|TT|)). The problem is that the fine-tuning grows 
very fast with M; in other words, the number of points 
in the parameter space with correct EW breaking de- 
cresases very quickly. In consequence this possibility is 
statistically penalized. On the contrary, for small M and 
large m, if tan /3 > 8, there is a focus-point region, with 
small fine-tuning. This region is statistically favoured, 



though this is counteracted by the penalization arising 
from the suppression in 5 SUSY a /i . This cannot be com- 
pensated by larger values of tan/3, since in this regime 
very big values of tan /3 (as would be needed for such 
compensation) start to be forbidden as we increase m. 
In consequence, a very large tan /3 is hardly favoured by 
an increasing m/,. 

Finally, let us mention that a lot of effort has been 
done in the literature to determine the most probable re- 
gion of the parameter space of the CMSSM [24 ] [27 1 [30 ] 
1321 1551 I59"H4"6] . This includes both Bayesian approaches 
(as the one followed here) and frequentist ones. The lat- 
ter (which can be considered as complementary to the 
Bayesian ones) are based on the analysis of the likeli- 
hood function in the parameter space. Thus they do 
not penalize regions from fine-tuning arguments (some- 
thing automatic in Bayesian analyses [27j[36] ). In con- 
sequence, following a frequentist approach it would be 
much more hard to show up the tension between rrih and 
g— 2. On the other hand, the present analysis differs from 
the previous ones in the fact that several hypothetic fu- 
ture scenarios, depending on the value of the Higgs mass, 
are considered and compared. 



V. CONCLUSIONS 

As it is well known, the SM prediction for the mag- 
netic anomaly of the muon, (basing the computation 
of the hadronic contribution on e + e~ data) shows a > 3 a 
discrepancy with the experimental result. It is common 
to consider this discrepancy as a signal of new physics 
(though, admittedly, the theoretical computation is con- 
troversial). In that case, SUSY is a most natural option 
for such new physics. 

However, as we have discussed in this paper, in the 
supersymmetric context there is a potential tension be- 
tween the requirement of SUSY contributions to the 
muon anomaly, <5 SUSY a AI , sufficient to reconcile theory 
and experiment, and a possibly large Higgs mass. In 
the CMSSM framework a large Higgs mass means O(10) 
GeV above the present experimental bound, rrih > 114.4 
GeV (for an SM-like Higgs). The tension arises because 
the main contributions to S MSSM a^ come from 1-loop 
diagrams with chargino-sncutrino and neutralino-smuon 
exchange, which grow with decreasing supersymmetric 
masses and increasing tan/3. But, on the other hand, 
in the MSSM the tree-level Higgs mass is bounded from 
above by Mz, so radiative corrections (which grow loga- 
rithmically with the stop masses) are needed to reconcile 
the theoretical predictions with the present experimen- 
tal lower bound. Thus, a large Higgs mass requires large 
supersymmetric masses, making impossible the task of 
reproducing the experimental value of a M . 

Although it is clear that in the limit of very large rrih 
(say above 135 GeV) the CMSSM must present the same 
discrepancy as the SM regarding the prediction for o^, 
it is much less clear for which size of rrih does the ten- 
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sion start to be unbearable, and would therefore put the 
model under pressure. Note that a particular value of the 
Higgs mass, say mh = 120 GeV, can be achieved through 
very different combinations of the supersymmetric pa- 
rameters, producing different values of <5 MSSM a A1 . On the 
other hand, it may happen that, for a given value of mh, 
the region of the parameter space compatible with a® xp is 
extremely tiny, implying a fine-tuning and thus a tension 
between the two observables. 

Our goal has been to quantify such tension, as a func- 
tion of mh, with the help of Bayesian techniques. As 
discussed at the end of sec. IV, this is the natural ap- 
proach if we want to incorporate the statistical penaliza- 
tion of fine-tuned regions of parameter-space. Certainly, 
if one just assumed a particular supersymmetric model 
(i.e. a point in the parameter-space, no matter how fine- 
tuned it were) then the statistical arguments used in this 
paper would not be appropriate. We have shown that 



for mh > 125 GeV the maximum level of discrepancy 
(~ 3.2 cr) is already achieved, indicating that SUSY has 
decoupled, and thus the prediction for a M coincides with 
the SM one. Given present-day data, requiring less than 
a 3 cr discrepancy, implies m/j < 120 GeV. This is a pre- 
diction of the CMSSM provided we accept the calculation 
of a M based on e + e~ data. For a larger Higgs mass we 
should give up either the CMSSM model (a the 3 cr level 
or above) or the computation of a M based on e + e _ ; or 
else accept living with such inconsistency. These are the 
main conclusions of this paper, and can be inferred di- 
rectly from Fig. 1. It is also important to note that, as 
discussed in section [TTH the CMSSM cannot remove the 
full 3.2cr discrepancy in a^. 

We have also examined the possibility that the ex- 
perimental uncertainty of ct° xp will decrease in the fu- 
ture, so that the discrepancy with the SM result be 5 cr, 
something that could happen in the next years. Then, 
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in the context of the CMSSM, a Higgs mass above 120 
GeV would imply a discrepancy larger than 4 a with 
the muon anomaly. Actually, the present lower bound, 
mil > 114.4 GeV, would already be inconsistent with the 
muon anomaly at the 3 a level. This illustrates the ten- 
sions within the CMSSM to accommodate a value of a p 
as the measured one (basing the theoretical calculation 
on present e + e~ data). 

Finally, we have shown how the probability distribu- 
tions of the most relevant parameters (universal scalar 
and gaugino masses, and tan fi) change with increasing 
77i h, which has obvious implications for the detection (or 
non-detection) of SUSY in the LHC. 
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